High T c superconductivity in superhard diamond-like BC5 
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Using density functional theory calculations we show that the recently synthesized [l| super- 
hard diamond-like BCs is superconducting with a critical temperature of the same order than that 
of MgB2. The average electron-phonon coupling is A = 0.89, the phonon-frequency logarithmic- 
average is (w)iog = 67.4 meV and the isotope coefficients are a(C) = 0.3 and a(B) — 0.2. In BC5, 
superconductivity is mostly sustained by vibrations of the B atom and its C neighbors. 

PACS numbers: 



Diamond is the hardest material available in nature, 
a very good thermal conductor and a large gap insula- 
tor. Heavily Boron-doped diamond (B^Ci-a) is super- 
conducting at Boron content 8 = 0.028 In thin films, 
T c increases with 5 up to 7K Q. These critical tem- 
peratures are lower than those of other superconducting 
carbon materials such as nanotubes i|| (T c = 15K), in- 
tercalated graphite compounds [1, H~_(T C = 11.4 K in 
CaC 6 ) and alkali-doped fullerides [J Q (T c = 33 K in 
RbCs2C6o). Thus, hole-doping of Diamond doesn't look 
very promising for high T c superconductivity. However, 
very recently, superhard nanocrystalline aggregates of 
BC5, corresponding to B-content of 5 — 0.166 were syn- 
thesized Besides being a superhard material, harder 
than c-BN, BC5 is interesting since the introduction of 
such a massive number of carriers in diamond can induce 
strong scattering between holes and hard phonon modes. 

At low B-content (S « 0.02 - 0.06), Density Func- 
tional TheoryfDFT) calculations based on a supercell 
approach (jj Till K)\ have shown that superconductivity 
is phonon mediated with holes scattering to high energy 
optical phonon modes. The contributions of the acoustic 
modes is negligible. Most surprisingly, B phonon states 
play a major role in the coupling to optical modes, despite 
the low B concentration. Moreover a Lowdin population 
analysis shows that the electronic states at e/ are mainly 
of B-character [9[. This is contrary to what happen in 
rigid band doping [l2| or in a Virtual Crystal Approach 



The coupling to hard optical phonons generates a large 
electron-phonon matrix element. However, at low dop- 
ing, the average electron-phonon coupling A and T c are 
small because of the low number of carriers. Indeed the 
increase of T c with doping observed in thin Films Q is 
mainly to due the increase in the density of states (DOS) 
[ici| at the Fermi level. 

These findings are valid at dopings of the order of some 
percents, but are questionable at dopings as large as those 
of BC5. First of all it is not clear at large doping how 
much the number of carriers can increase, since a rigid 
band picture fails and the results depends on the fine de- 
tails of the DOS. Furthermore, if many electronic states 



are allowed to couple to phonons generating a consider- 
able electron-phonon coupling, the system can be driven 
to a charge density wave insulator. Thus the occurrence 
of superconductivity in BC5 cannot be inferred from low 
doping results. 

In this work, using density functional theory calcula- 
tions, we demonstrate that BC5 has a critical tempera- 
ture of the same order of that of MgB2 . 

Cubic BC5 occurs as nanocrystalline aggregates 
The crystal structure has, on average, cubic symmetry 
with volume per atom of 6.00A 3 . This value is 6% larger 
than the diamond volume per atom, 5.67A 3 . Due to their 
similar atomic numbers, diffraction cannot distinguish 
between B and C atoms and consequently the position 
of B in the cell is not determined. In order to resolve the 
BC5 crystal structure we start with the 6-atoms hexag- 
onal supercell of the 2 atoms diamond cell. The hexag- 
onal supercell is compatible with a cubic symmetry for 
c/a = \f& as 2.45. The diamond theoretical equilibrium 
structure has a = 2.529A Then we replace a C atom 
with a B atom and we perform volume and force opti- 
mization [lij]. The results are reported in tab. Q] We 
find a = 2.55 A and c/a = 2.50. The theoretical vol- 
ume per atom in BC5 is 5.2% larger than in diamond, in 
very good agreement with experiments. The c/a found is 
only slightly larger than the ideal one for cubic symmetry 
so that the BC5 cell can be seen as a small elongation of 
the cubic cell along the cubic (111) axis. The most stable 
structure is non-magnetic. Then we consider a 12 atoms 
cell obtained replicating the hexagonal cell along c and 
consider all the possible positions for two B atoms in this 
supercell. We found the configuration corresponding to 
the 6 atoms hexagonal unit cell to be the most stable by 
5 mRyd/atom. 

The BC5 and diamond (Ce) electronic structures along 
high energy directions in the hexagonal Brillouin zone 
are plotted plotted in Fig. [TJ For Cq we used the BC5 
lattice parameters and the Fermi level refers to a rigid 
band doping of diamond of 1 hole per Cq unit. As it 
can be seen, the electronic structure of BC5 cannot be 
interpreted as rigid band doping of diamond. 

In BC5, the Fermi level e/ is shifted of 2.44 eV respect 
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Atom type 



X 



Y 



B 

C 
C 

c 
c 
c 



0.0 
0.0 

1/3 
1/3 
2/3 
2/3 



0.0 
0.0 

-1/3 
-1/3 
-2/3 
-2/3 



0.0070 
0.2592 
0.3386 
0.5816 
0.6646 
0.8999 



TABLE I: Reduced coordinates of the theoretically devised 
atomic structure of BC5 . The crystal symmetry is hexagonal 
with a = 2.55A and c/a = 2.50. 



density wave instabilities occur in the system. Decom- 
position of phonon vibrations (Fig. [2] ) into atomic 
components shows that, despite B being lighter than C, 
the harder phonon modes (120-150 meV) are due to C- 
vibrations. At lower energies (< 100 meV) the B com- 
ponent in the phonon density of states (PHDOS) is from 
1/3 to 1/2 of the C-component. 

The superconducting properties can be understood cal- 
culating the electron-phonon coupling A ql/ for a phonon 
mode v with momentum q and phonon frequency w q „, 
namely: 



to the top of the valence band, to be compared with a 
shift w 0.8 eV at 6 = 0.0278 0. The DOS at e/ is 
N(0) = 0.78 states/eV/(6 atoms cell), 2.1 times larger 
than its value at S — 0.0278 [nj, meaning that the num- 
ber of carriers induced by B-doping grows continuously 
without saturating even for such large S. Lowdin popu- 
lation analysis demonstrates that at e/ the total B DOS 
is Nb(0) = 0.19, states/eV/(6 atoms cell), slightly larger 
that N(Q)/6. Thus, contrary to what happens at low 
doping [9J] , the number of carriers in B electronic states 
is comparable to the B-content. 
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FIG. 1: BC5 and diamond (Ce) band structure and DOS in 
the hexagonal Brillouin zone using the BC5 lattice parame- 
ters. The DOS is in units of states/eV/(6 atoms cell). In Ce 
the Fermi level corresponds to rigid band doping of diamond. 
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the sum is performed over a 
mesh in the Brillouin Zone. The 

5kn,k+qm = (kn|<SV/<Je q „|k + qm)/^, 
is the Kohn-Sham potential and e q , y = 
q ' ■ ' is the Fourier transform 



is 
V 



of the a component of the phonon displacement of the 
atom A in the unit cell, Ma is the mass of atom A and 
e\ are Aa components of phonon eigenvector nor- 
malized on the unit cell. 

The average electron-phonon coupling is A = 
S ql y ^qy/N q ~ 0.89. Thus BC5 is a moderate coupling 
superconductor with A comparable to MgB 2 . To ad- 
dress the role of the different vibrations in determining 
the electron-phonon coupling we decompose the electron- 
phonon coupling [l^ | into selected atomic vibrations, us- 
ing the relation: 



A 



A,B 



q\Aa,B0 



[Cc 



\Bf3,Aa 



a,p * q 

(2) 

where A, a (Bf3) indicates the displacement of the 
A th (B th ) atom in a (/3) Cartesian direction, [GJao.b^ = 
J2^ m 4 9* Aa 9Bf38(e kn )S(e k+(lm )/[N(0)N k }, and g Aa = 

(kn\SV/5uqAa\k + qm)/v2- The C q matrix is the 
Fourier transform of the force constant matrix (the 
derivative of the forces with respect to the atomic dis- 
placements). The decomposition [20[ leads to the follow- 
ing matrix (A = J2i , 




B C C" 

0.50 -0.18 0.04 

-0.18 0.69 -0.07 

0.04 -0.07 0.12 



(3) 



The BC5 phonon dispersion [l6| is shown in Fig.O No 
dynamical instabilities are seen in the phonon spectrum 
meaning that the crystal structure obtained with geo- 
metrical optimization is dynamically stable. No charge 



where C stands for the first two C nearest neighbors of 
the B atom, and C" stands for the other three C atoms 
in the cell. The dominant contribution to the electron 
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FIG. 2: (Color). (Left) BC5 phonon dispersion in the hexagonal Brillouin zone. The radius of the red dots is proportional to 
the B content of the branch. (Right) BC5 Phonon density of states projected over B and C states, Eliashberg function a 2 F(u>) 
and integrated electron-phonon coupling A(w). The Eliashberg function has been decomposed three different contributions, due 
to vibrations of B and its C neighbors (BC), and vibrations of far neighbors (C"). 



phonon-coupling comes from the vibrations of the B atom 
and its C neighbors. 

The Eliashberg function 



q'' 



(4) 



and the integral A(w) = 2 f£ dcu'a 2 F(uj')/uj' are shown 
in Fig. [2] . a 2 F(uj) is composed by three main peaks. 
This is different from the low doping case ll| since in 
BC5 the high energy modes (80-90 meV and 130 meV) 
are softer due to higher doping and the low energy modes 
(« 40 meV) are significantly coupled (at low doping the 
coupling is very weak). 

To understand the role of different vibrations in deter- 
mining the total electron-phonon coupling, we decompose 
a 2 F as: 



a 2 F{u) = J2 a XsF(Lo) 
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(e£^)*. Then we plot 



where [L q ] B /3,Aa = E P e qP Wqp v- qP 
the decomposition of a 2 F(u>) into vibrations of B and its 
first two C nearest neighbors (BC) and vibrations of the 
other three C atoms in the cell (C"). As it can be seen 



in Fig. [21 the contribution of the highest energy mode 
(115-160 meV), essentially due to vibrations of the C-C 
bonds, to the electron-phonon coupling is 1=3 0.13. This is 
half of the contribution of the high energy structure to A 
in doped-diamond at S = 0.018 (from the inset of fig. 3 in 
ref. [ll[ one can infer that the contribution of the lower 
energy modes is « 0.09), indicating that in BC5 substan- 
tial weight is shifted at low energy. At lower energies, 
the other two main peaks in a 2 F(oj) are entirely due to 
B and its C neighbors. Further decomposition of a 2 F(ui) 
onto vibration of B and C atoms (not shown) indicates 
that the two low energy peaks are due to phonon modes 
with concerted movements of BC atoms. Thus, despite 
the lowest B concentration, the contribution of vibrations 
associated to the B atom and its neighbors are the main 
responsible for superconductivity. 

As demonstrated, BC5 has two times larger number 
of carriers than doped-diamond at S = 0.0278 and the 
phonons coupled to electrons are less energetic. Since A ql/ 
is proportional iV(0)/w q;/ , both these effects cooperates 
in increasing A. 

The critical superconducting temperature is estimated 
using the McMillan formula [2l|, 



T c = 



(w)log 



exp 



1.2 

where fj,* is the screened ^uuiumu pt 
(w)iog = exp[f J + °° a 2 F(cj) log(w)/w 



(6) 



1.04(1 + A) 
A-/i*(l + 0.62A) J 

Coulomb pseudopotential and 
,.,\ ir^a-d.A l>.>du)] is the phonon 



4 



frequencies logarithmic average. We obtain (w)i og = 67.4 
meV, to be compared with 105 meV in doped diamond at 
6 = 0.027810] and 62 meV in MgB 2 [22j]. The reduction 
in (cj)i g as compared to hole doped diamond is due to (i) 
activation of coupling to low energy modes which are very 
weakly coupled in B-doped diamond and (ii) softening of 
the high energy modes due to the larger doping. 

Using the same value of fi* = 0.1 neces sary for hole- 
doped diamond in order to obtain T c = 4K [l0( we obtain 
T c = 45 K, which is larger than the experimentally mea- 
sured 39 K in MgB 2 [15| and puts BC5 in the class of 
high T c superconductors. As in MgB 2 the T c could be 
further increased by multi-band effects [2 31] 

The mutual relevance of B and C phonon modes in sus- 
taining superconductivity in BC5 can be addressed mea- 
suring the isotope effect coefficients for a given atomic 
specie X = B,C, namely a(X) = -(dlogT c /dM x ). We 
obtain a(C) = 0.3 and a(B) = 0.2, confirming the im- 
portant role of B-phonon modes. 

A question arise whether the actual samples of BC5 are 
able to sustain superconductivity due to the reduced size 
of the grains (diameter 10-15 nm) [l[. The observation of 
Bulk superconductivity is possible only if if the coherence 
length £0 ~ hvp/A is at least com par able to the size 
of the grains. In B-doped diamond [24|, £0 = 240 nm. 
Assuming parabolic bands and the same A(0)/kt,T c = 
1.78 as in Boron-doped diamond, v/ is 16 times larger and 
£0 ~ 400 nm. Thus it is necessary to grow larger samples 
to observe bulk superconductivity in BC5, possibly by 
longer synthesis or by the use of catalysts to speed up 
the reaction. 

High energy phonon modes can lead to large T c s even 
with moderate electron-phonon scattering (A ~ 0.8 — 1.0), 
as in MgB 2 . Interestingly, the quest for such a kind of 
high-temperature phonon-mediated superconductors co- 
incides with that of metallic superhard materials. Su- 
perhard materials have large elastic constants, requiring 
short and strong chemical bonds, typically found in light- 
element compounds. These conditions result in energetic 
phonon modes. Thus high T c superconductivity could be 
realized in superhard materials if doping is large enough 
to sustain a moderate electron-phonon coupling A. 

In this work we have shown that this happens in 
the recently-synthesized 4j superhard BC5 which is pre- 
dicted to be metallic and superconducting with a T c = 45 
K, the largest T c ever for a phonon-mediated supercon- 
ductor. 
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